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ABSTRACT:  The  racemic  product  of  the  Betti  reaction  of 
5-chloro-8-hydroxyquinoline^  benzaldehyde^  and  2-aminopyri- 
dine  was  separated  by  chiral  HPLC  to  determine  which 
enantiomer  inhibited  botulinum  neurotoxin  serotype  A.  When 
the  enantiomers  unexpectedly  proved  to  have  comparable 
activity^  the  absolute  structures  of  (+)-(R)-1  and  (— )-(S)-l 
were  determined  by  comparison  of  calculated  and  observed 
circular  dichroism  spectra.  Molecular  modeling  studies  were 
undertaken  in  an  effort  to  understand  the  observed  bioactivity 
and  revealed  different  ensembles  of  binding  modes^  with 
roughly  equal  binding  energies^  for  the  two  enantiomers. 


Cl 


KEYWORDS:  Chiral  resolution^  Betti  reaction  products^  TDDFT  CD  calculations^  molecular  dockings  inhibition  of  botulinum 
neurotoxin 


We  recently  reported  the  identification  of  7-substituted 
8-hydroxyquinolines  (e.g.^  l);  products  of  the  Betti  reac- 
tiou;  as  leads  to  potential  chemotherapies  for  botulinum  poison¬ 
ing.^  The  Betti  reaction^  originated  as  a  condensation  of  Schiff 
bases  with  2-naphthol;^'^  but  it  was  later  extended  to  other 
nucleophilic  aromatic  substrates^  e.g._,  8-hydroxyquinoline.^ 
While  a  chiral  center  is  formed  in  this  reaction^  the  products 
are  typically  racemic  mixtures^  as  the  reagents  are  usually  achiral 
(Scheme  l).  Over  the  years^  there  have  been  several  efforts  to 
resolve  Betti  product  racemates^  most  of  which  involve  preparing 
diastereomeric  salts  with  chiral  acids.^”^ 

Chiral  chromatography  seemed  an  obvious  alternative  for 
resolution  of  these  racemates^  but  we  examined  a  total  of  six 
chiral  columns  and  18  different  methods  before  observing 
sufficient  resolution  of  the  two  enantiomers  on  a  Chiralcel  OD 
column  to  permit  semipreparative  purification  of  adequate 
quantities  of  (+)-l  and  (— )-l  for  evaluation  in  our  bioassays 
and  assignment  of  absolute  configuration  (Figures  SI  and  S2  of 
the  Supporting  Information). 

Initial  comparison  of  the  botulinum  neurotoxin  serotype  A 
(BoNT/A)  inhibitory  activity  of  the  (+)  and  (— )  enantiomers  of 
1  was  accomplished  via  an  HPLC-based  assay  using  a  full-length 


recombinant  BoNT/A  protease  light  chain  (LC).^  To  our 
surprise^  both  compounds  displayed  similar  IC50  values^  I.O 
and  1. 1  juM  for  (+)-l  and  (—)-!;  respectively.  We  subsequently 
evaluated  their  potential  to  inhibit  the  biological  activity  of 
BoNT/A  holotoxin  in  murine  neuroblastoma  N2a  cells. ^  No 
difference  was  observed  in  percent  inhibition  (P  >  0.05)  of 
SNAP-25  cleavage  for  both  enantiomers  and  the  racemate  (±)-l 
at  the  four  concentrations  tested  {60,  AS,  30^  and  15  fiM).  We 
then  examined  the  efficacy  of  these  compounds  in  mouse  phrenic 
nerve  hemidiaphragm  preparations  (MPNHDA).^  Similar  to 
observations  in  HPLC  and  cell-based  assays^  both  (+)-l  and 
(— )-l  were  equipotent  (P  =  0.94)  in  the  tissue-based  assay.  At  2 
fiM  concentrations^  both  enantiomers  dramatically  delayed  (P  = 
1.58  X  I0~^  and  2.30  x  I0~^  for  (+)-l  and  (— )-l;  respectively) 
the  BoNT/A-induced  paralytic  half-time  3-fold.  The  compara¬ 
tive  testing  of  (±)-l  and  the  two  enantiomers  is  summarized  in 
Table  I. 
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We  then  turned  to  assigning  the  absolute  configuration  of 
(+)-l  and  (— )-l  via  comparison  of  calculated  and  experimental 
electronic  dichroism  (CD)  spectra.  As  a  prelude^  we  determined 
the  3-dimensional  conformation  of  1  through  a  series  of  NMR 
experiments.  Proton  and  carbon  resonances  were  assigned  from  a 
combination  of  COSY_,  HSQC;  and  HMBC  experiments  (Table 
Sh  Supporting  Information).  Numerous  NOE  interactions  were 
observed  in  the  NOESY  experiment  (Scheme  2),  with  those  of 
NH;  OH;  H6;  and  H9  providing  the  greatest  insight  into  the 
conformation  of  the  compound.  While  all  NMR  data  were 
obtained  for  the  racemic  mixture;  the  (S) -enantiomer  is  shown 
for  ease  of  portraying  the  structural  analysis.  The  observed  NOE 
between  H6  and  H18a;b  and  lack  of  NOE  between  the  C8 
hydroxyl  proton  and  H18a;b  indicate  that  the  phenyl  ring  is 
rotated  away  from  the  C8  OH.  The  NOE  between  the  OH  and 
the  NH;  as  well  as  the  OH  and  H9  interaction;  supports  this 
conformation.  The  phenyl  ring  assumes  a  position  perpendicular 
to  the  quinoline  ring  system  to  minimize  steric  interactions.  This 
is  supported  by  NOEs  of  H18a;b  with  H6;  H9;  and  the  NH. 
Interestingly;  HIS  and  HI6  both  have  an  NOE  with  the  NH;  but 
not  with  H9;  suggesting  that  the  pyridine  ring  is  oriented  so  as  to 
place  the  ring  nitrogen  toward  H9;  while  HIS  and  HI 6  are 
projected  away.  AlsO;  this  allows  the  pyridine  ring  to  be  ortho¬ 
gonal  to  the  quinoline  and  pseudoparallel  to  the  phenyl  ring; 
further  minimizing  steric  interactions.  A  locked  conformation  of 
the  compound  due  to  hydrogen  bonding  between  the  lone  pair  of 
the  NH  and  the  proton  of  the  C8  hydroxyl  is  highly  probable. 
However;  the  high  number  of  NOEs  for  the  NH  (OH;  H6;  H9; 
HI8a;b;  and  HI9a;b)  implies  that  the  structure  may  not  be 
entirely  rigid. 

To  assign  the  absolute  configuration  of  (+)-  and  (— )-l^  their 
electronic  CD  spectra^'^^  were  recorded  in  solution  and  com¬ 
pared  with  those  calculated  using  the  time-dependent  DFT 


Scheme  1 


(TDDFT)  method. The  enantiomers  of  1  show  almost 
mirror  image  CD  spectra  in  methanol;  as  expected  (Figure  l). 
Due  to  the  presence  of  three  different  aromatic  chromophoreS; 
these  spectra  feature  many  bands  between  200  and  3S0  nm.  For 
the  (+) -enantiomer;  the  first  band  appearing  in  the  280—340  nm 
region  is  broad  and  positive;  followed  by  a  moderately  intense 
negative  band  centered  around  260  nm  and  two  stronger  bands;  a 
positive  one  at  245  nm  and  a  negative  one  around  220  nm. 

Several  rotatable  bonds  in  the  structure  of  1;  most  directly 
affecting  the  relative  orientation  of  the  chromophoreS;  made 
obtaining  a  reliable  set  of  input  structures  crucial  for  CD 
calculations.  A  preliminary  molecular-mechanics  conformational 
search  was  run  with  the  MMFF  force  field;  using  a  starting 
geometry  with  (R)  absolute  configuration  (see  Supporting 
Information  for  details).  All  low-energy  structures  obtained  were 
optimized  with  the  DFT  method  at  the  B3LYP/6-3IG(d)  level; 
converging  to  a  set  of  nine  distinct  conformers  within  10  kj/mol. 
Their  energies  and  populations  at  298.15  K  were  estimated  with 
B3LYP/6-3II++G(d;p)  in  methanol  (PCM  solvent  model). 
The  low-energy  DFT  structures  were  then  checked  against  the 
conformational  picture  provided  by  NMR  experiments;  bearing 
in  mind  that  the  coexistence  of  several  low-energy  minima 
rendered  interpretation  of  NOEs  in  terms  of  a  single  conforma¬ 
tion  questionable.  The  calculated  structures  may  be  divided  into 


Scheme  2 


Table  1.  Comparative  Testing  of  (=t)-l;  (+)-l,  and  (—)-!'* 

HPLC  assay  N2a  cell  culture  assay  MPNHDA^ 


%  inhibition  IC50  %  inhibition  of  SN7\P-25  cleavage  minutes^^/ 


sample 

20  piM 

5  piM 

CaM) 

60  piM 

45  piM 

30 /iM 

15  piM 

2  piM 

(±)-i 

92 

88 

1.5 

91 

90 

74 

41 

ND" 

(+)-i 

94 

90 

1.1 

96 

94 

89 

56 

I9E 

(-)-i 

91 

90 

1.0 

96 

93 

69 

54 

188^ 

^Assays  conducted  as  described  in  ref  1  and  the  Supporting  Information.  ^ Mouse  phrenic  nerve  hemidiaphragm  assay.  ‘^Average  time  to  50%  loss  of 
twitch  tension  (min).  ^Average  value  for  BoNT/A  toxin  control  was  63  min.  ^  Not  determined. value  <  0.001  (highly  significant)  for  comparison  with 
values  recorded  for  the  BoNT/A  control;  statistical  analyses  performed  using  SigmaPlot  10  (Systat  Software;  San  Jose;  CA). 
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two  subsets  according  to  the  rotation  around  the  C7— C9  bond. 
In  one  (major)  subset^  composed  of  four  lowest-energy  con- 
formers;  H9  is  directed  toward  the  OH.  Within  this  subset^  two 
conformers  are  especially  stable  (absolute  minimum^  30.8% 
population  at  298.15  K;  second  lowest  minimum^  +0.37  kj/ 
mol;  26.5%;  Figure  2);  the  structures  and  relative  energies  for  all 
minima  are  reported  in  the  Supporting  Information.  In  the 
second  (minor)  subset;  composed  of  the  remaining  five  con- 
formers  accounting  for  25%  overall  population;  H9  is  again  in  the 
plane  of  the  quinoline  ring  but  directed  toward  C6.  These 
findings  corroborate  the  observed  NOEs  that  NH  and  H9  each 
have  with  both  the  C8  OH  and  H6. 

Taking  the  lowest-energy  structure  as  a  test  molecule; 
TDDFT  calculations  were  run  using  different  combinations  of 
hybrid  DFT  fimctionals  (B3LYP;  PBEO;  CAM-B3LYP;  BH& 
HLYP)  and  basis  sets  (SVP;  TZVP;  aug-TZVP);^"^  in  vacuo  or  in 
methanol.  The  three  basis  sets  led  to  very  similar  results  (using 
B3LYP);  except  for  a  small  red  shift  observed  for  all  computed 
transitions  on  increasing  the  basis  size.  It  appears  that  the  two 
“standard”  functionals  B3LYP  and  PBEO  led  to  calculated 
transition  energies  underestimated  with  respect  to  the  experi¬ 
ment;  while  the  opposite  was  true  for  Coulomb-attenuated 
B3LYP  (CAM-B3LYP)  and  the  half  a  d-half  functional 
BH&HLYP.  Looking;  for  example;  at  the  absorption  band 
measured  at  250  nm;  the  calculated  transition  wavelength  was 
270  and  230  nm  with  B3LYP/aug-TZVP  and  CAM-B3LYP/ 
SVP;  respectively.  Apart  from  a  systematic  wavelength  shift;  the 


200  250  300  350 


A,/  nm 

Figure  1.  CD  spectra  of  the  enantiomers  of  1  in  methanol;  concentra- 


shape  of  the  calculated  CD  spectrum  was  similar  in  all  cases. 
Finally;  including  the  solvent  model  in  CAM-B3LYP/SVP 
calculations  did  not  appreciably  change  calculated  frequencies 
and  spectra. 

TDDFT  calculations  were  then  run  on  all  low-energy  struc¬ 
tures  using  the  two  combinations  B3LYP/ aug-TZVP  and  CAM- 
B3LYP/SVP  in  vacuo.  The  resulting  CD  spectra  were  weighted 
with  the  respective  Boltzmann  factors  (estimated  from  B3LYP/ 
6-3II+G(d;p)  internal  energies  in  methanol)  at  298.15  K 
and  summed  to  afford  weighted  average  spectra.  In  all  caseS; 
input  structures  had  (R)  absolute  configuration.  Figure  3  displays 
spectra  computed  with  the  two  methods  for  the  lowest-energy 
structure  and  the  weighted  averages  over  nine  structures.  Apart 
from  the  already  discussed  wavelength  shift  (taken  into  account 
by  a  frequency  correction  in  Figure  3);  the  overall  shapes  of  the 
four  spectra  are  quite  similar;  especially  in  the  low-energy  region 
(where  TDDFT  calculations  are  intrinsically  more  accurate). 
CAM-B3LYP/SVP  results  agree  especially  well  with  the  experi¬ 
mental  spectrum  for  (+)-l  in  sign;  position;  and  intensity  of 
bands  (Figure  3b).  Therefore;  the  absolute  configuration  of  the 
enantiomers  of  compound  I  maybe  assigned  as  (+)-(R)-1  and 
(— )-(S)-l.  It  must  be  noted  that  the  calculated  average  CD 
spectrum  is  the  superposition  of  very  heterogeneous  component 
spectra;  thuS;  the  apparent  bands  are  actually  the  convolution 
from  several  transitions  (Supporting  Information). 

The  unexpected  observation  of  virtually  equivalent  BoNT/ 
A-inhibitory  activity  for  the  enantiomers  in  three  different 
bioassays  prompted  us  to  examine  the  bound  state  of  the  ligands 
to  rationalize  the  apparent  lack  of  discrimination.  Given  the 
observed  pharmacological  activity;  we  hypothesized  that  both 
enantiomers  bind  BoNT/A;  presumably  in  the  Zn^^-containing 
active  site.  Nevertheless;  experimental  observations  suggest  that 
binding  of  a  similar  ligand  is  not  exclusively  due  to  Zn^^ 
chelation. To  test  this  hypothesis  computationally;  we  as¬ 
sessed  relative  binding  free  energies  for  the  enantiomers  by 
performing  Thermodynamic  Integration  with  Molecular  Dy¬ 
namics  simulations^^  (see  the  Supporting  Information  for 
details).  Because  the  protonation  state  of  the  bound  ligand  is 
unknown;  we  assumed  that  both  enantiomers  bind  the  active  site 
Zn^^  without  ionization  of  the  OH  group.  Using  10  lowest  stable 
conformations  of  each  ligand;  we  generated  corresponding 
complex  conformations  by  superposing  the  quinoline  moiety 
onto  a  template  complex  derived  with  AutoDock4.0.^^~^^  The 
resulting  calculations  indicated  that  the  (—)-S-enantiomer  would 


Figure  2.  DFT-optimized  structures  for  the  absolute  lowest  energy  (left)  and  second  lowest  energy  (right)  conformers  of  (-R)-l. 
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Figure  3.  CD  spectra  calculated  for  (-R)-l  with  two  TDDFT  methods:  (A)  B3LYP/aug-TZVP;  (B)  CAM-B3LYP/SVP.  Blue  lines:  spectra  calculated 
for  the  lowest-energy  DFT  structure  (divided  by  2  for  better  comparison).  Black  lines:  averages  of  spectra  calculated  for  nine  low-energy  DFT  structures 
weighted  with  Boltzmann  factors  at  298.15  K  using  populations  estimated  with  B3LYP/6-311G++(d,p)  in  methanol.  Red  dashed  line:  experimental 
spectrum  for  (+)-l.  Frequency  corrections  of  +2500  cm~^  and  —2000  cm~^  have  been  applied  to  calculated  spectra  in  panels  A  and  B^  respectively. 


A 


Figure  4.  Examples  of  binding  modes  for  the  (R) -enantiomer  (A)  and 
(S) -enantiomer  (B).  Active  site  residues  (yellow)  of  the  protein  (green) 
participate  in  binding  of  the  Zn^^  (magenta)  and  the  ligand  (gray).  Blue 
and  red  colors  correspond  to  oxygen  and  nitrogen  atoms.  Certain  active 
site  residues  and  hydrogen  atoms  are  omitted  for  clarity.  The  red  dashed 
lines  show  some  coordinating  and  hydrogen  bonds. 


bind  BoNT/A  by  3.8  zb  2.4  kcal/mol  more  favorably  than  its 
( + )  -R-counterpart. 

The  large  standard  deviation  in  the  free  energy  calculations 
indicates  significant  variability  in  the  conformational  ensembles 
of  the  bound  ligands.  To  analyze  the  variety  of  binding  modes  of 


both  (+)-R-l  and  (  — )-S-l  to  BoNT/A^  along  vfith  the  corre¬ 
sponding  binding  features^  we  performed  GROMOS-style 
clusterings^  of  the  sparsely  saved  simulation  trajectories  in  the 
physical  state  for  each  of  the  10  starting  configurations  of  the  (R)- 
and  (S)-  complexes.  (Table  S3  of  the  Supporting  Information 
summarizes  the  16  most  prominent  features  describing  the  two 
ensembles  and  comprising  H-bonding  interactions^  coordinating 
bonds  to  the  Zn^^  ion  and  7t—7t  interactions.)  The  most  robust 
features  include  coordinate  bonds  of  atoms  N1  (quinoline)^ 
and  N12  (pyridine)  to  the  Zn^^  ion.  Depending  on  the  number 
of  bonds;  the  ligand  can  be  monodentate  or  bidentate.  When 
both  O  and  N12  chelate  the  Zn^^_,  the  NIO  (amine)  atom  is 
brought  sufficiently  close  to  the  Zn^^  that  this  coordination 
mode  could  be  considered  as  tridentate.  The  most  common 
H-bonding  interaction  is  that  of  the  coordinated  OH  group  with 
the  carboxyl  of  E224. 

Clustering  analysis  of  the  trajectories  of  the  protein— ligand 
complexes  captured  the  diversity  of  the  binding  modes  that  can 
be  described  by  a  few  features.  Most  of  these  features  are 
common  to  both  enantiomers  and  pertain  to  specific  interactions 
of  the  ligand  with  the  binding  site  residues.  We  found  that  the 
(— )-S-l  ensemble^  on  average^  has  more  interaction  features  per 
cluster^  3.6  versus  3.0  for  that  of  the  (+)-R-l.  In  addition^  the 
(— )-S-l  ensemble  has  a  larger  number  of  unique  interactions 
with  the  active  site  compared  to  that  of  the  (+)-R-l.  These 
observations  are  consistent  with  the  computed  binding  free 
energy  difference.  Figure  4  illustrates  examples  of  the  likely 
binding  mode  for  each  enantiomer. 

Only  a  few  interactions  exist  that  can  differentiate  between  the 
two  enantiomers.  In  particular^  the  NlOH— E224  interaction  is 
realized  only  by  the  (R)-enantiomer;  whereas  the  interactions 
OH-Y366;  N12— F163(NH);  and  Ar:Phe-Y366  are  only  realized 
by  the  (S) -enantiomer.  Although  it  is  possible  that  different 
realizations  of  common  interaction  features  by  the  enantiomers 
can  further  contribute  to  binding  disparity^  quantifying  those 
differences  could  be  daunting. 

The  fact  that  our  calculations  favor  BoNT/A  binding  by  the 
(— )-S-l  suggests  that;  in  reality^  the  ligands  might  bind  differ¬ 
ently.  Unlike  experimental  measurements^  our  calculated  values 
provide  a  more  limited  assessment  of  the  relative  binding 
affinities  of  the  two  enantiomers  under  specified  assumptions. 
It  is  also  possible  that  the  ligand  can  undergo  epimerization; 
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similar  to  that  of  ibuprofen^^'^^  in  the  presence  of  certain 
enzymes. 

In  conclusion^  we  separated  the  enantiomers  of  1  (chiral 
HPLC)^  established  their  chemical  structure  (NMR)  and  abso¬ 
lute  configuration  (CD;  molecular  modeling);  evaluated  their 
BoNT/A  inhibitory  activity;  and  explored  their  docking  motifs 
with  BoNT/A  LC.  To  our  knowledge;  this  is  the  first  use  of  CD 
calculations  to  assign  the  absolute  configuration  of  Betti  reac¬ 
tion  product  enantiomers.  While  chiral  chromatography  has 
recently  been  applied  to  the  separation  of  naphthol-based  Betti 
productS;^^'^^  this  is  the  first  such  separation  of  Betti  products 
comprised  of  8-hydroxyquinoline;  an  aryl  aldehyde  and  an 
aryl  amine. 

While;  in  a  vast  majority  of  caseS;  one  enantiomer  of  a  racemic 
drug  or  drug  candidate  has  significantly  greater  pharmacological 
activity  than  the  other;  as  has  recently  been  demonstrated  for  one 
BoNT/A  inhibitor;^^  we  found  essentially  equivalent  BoNT/ 
A-inhibitory  activity  in  (+)-(R)-1  and  (— )-(S)-l.  However;  this 
unexpected  finding  can  be  explained  by  the  proposed  docking 
motifs  for  the  two  enantiomers — different  ensembles;  but  nearly 
equivalent  in  energy.  We  are  currently  seeking  to  confirm  those 
binding  models  by  crystallizing  each  enantiomer  in  the  active  site 
ofBoNT/A. 
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